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Abstract 

We compare 'fixed flavor number scheme' (FFNS) and 'variable flavor number 
scheme' (VFNS) parton model predictions at high energy colliders. Based on our 
recent LO- and NLO-FFNS dynamical parton distributions, we generate radiatively 
two sets of VFNS parton distributions where also the heavy quark flavors h = c,b,t 
are considered as massless partons within the nucleon. By studying the role of these 
distributions in the production of heavy particles (hh, tb, hW ± , Higgs-bosons, etc.) 
at high energy ep, pp and pp colliders, we show that the VFNS predictions are com- 
patible with the FFNS ones (to within about 10-20% at LHC, depending on the 
process) when the invariant mass of the produced system far exceeds the mass of the 
participating heavy quark flavor. 



In a recent publication [1\ we updated the dynamical leading order (LO) and next-to- 
leading order (NLO) parton distributions of [2J. These analyses were undertaken within the 
framework of the so called 'fixed flavor number scheme' (FFNS) where, besides the gluon, 
only the light quarks q = u,d,s are considered as genuine, i.e. massless partons within 
the nucleon. This factorization scheme is fully predictive in the heavy quark h = c,b,t 
sector where the heavy quark flavors are produced entirely perturbatively from the initial 
light quarks and gluons - as required experimentally, in particular in the threshold region. 
However, even for very large values of Q 2 , Q 2 ^> m 2 b , these FFNS predictions are in 
remarkable agreement with DIS data [HE] and, moreover, are perturbatively stable, despite 
the common belief that 'non-collinear' logarithms ln(Q 2 /m 2 ) have to be resummed. 

In many situations, however, calculations within this factorization scheme become un- 
duly complicated. For example, the single top production process at hadron colliders 
via VT-gluon fusion requires the calculation of the subprocess ug — > dtb at LO and of 
ug — ► dtbg etc. at NLO. It thus becomes expedient to consider for such calculations the 
so called 'variable flavor number scheme' (VFNS) where also the heavy quarks h — c,b,t 
are taken to be massless partons within the nucleon. In this scheme, the above FFNS 
calculations simplify considerably, i.e. one needs merely ub — > dt at LO and ub — > dig 
etc. at the NLO of perturbation theory [3]. The VFNS is characterized by increasing 
the number n/ of massless partons by one unit at Q 2 = m\ starting from n/ = 3 at 
Q 2 = m 2 , i.e. c{x, m 2 ) = c(x,m 2 ) = 0. The matching conditions at LO and NLO are 
fixed by continuity relations [4j at the respective thresholds Q 2 = m\. Thus the 'heavy' 
rif > 3 quark distributions are perturbatively uniquely generated from the Uf — 1 ones via 
the massless renormalization group Q 2 -evolutions. The running strong coupling can be 
approximated by the common NLO 'asymptotic' solution 

a s (Q 2 ) 1 ft lnln(Q 2 /A 2 ) 

4tt " /? ln(Q 2 /A 2 ) (31 [ln(Q 2 /A 2 )] 2 U 

with p = 11 — 2n//3 and f3\ = 102 — 38n//3, which turns out to be sufficiently accurate pQ 
for our relevant Q 2 -region, Q 2 ,> 2 GeV 2 . Since /?o,i are not continuous for different flavor 
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numbers n/, the continuity of a s (Q 2 ) requires to choose different values for the integration 
constant A for different rif, A^ n f\ which are fixed by the common a s (Q 2 ) matchings at the 
flavor thresholds Q = m^t- Choosing m c = 1.3 GeV, mj = 4.2 GeV and m t = 175 GeV, 
one obtains A^' 6) = 269.7, 184.5, 72.9 MeV according to our dynamical NLO(MS) fit [lj 
which resulted in a s [M 2 ) = 0.1145. In LO, where ft = 0, we obtained [TJ A^ 5 ' 6) = 181.8, 
138.3, 70.1 MeV corresponding to a s (M|) = 0.1263. 

Our choice for the input of the 'heavy' VFNS distributions are the LO and NLO dy- 
namical FFNS distributions [lj at Q 2 = m 2 . The VFNS predictions at scales Q 2 ^> m 2 
should become insensitive to this, somewhat arbitrary, input selection [5] whose virtues are 
the fulfillment of the standard sum-rule constraints together with reasonable shapes and 
sizes of the various input distributions. As we shall see, this expectation is based on the 
fact that at Q 2 ^> ml the VFNS distributions are dominated by their radiative evolution 
rather than by the specific input at Q 2 = m\. In other words, because of the long evolution 
distance, input differences get 'evolved away' at Q 2 ^> m\ where the universal perturbative 
QCD splittings dominate. 

As a first test of the VFNS 'heavy' quark distributions we consider charm and bottom 

electroproduction processes, since deep inelastic structure functions play an instrumental 

role in determining parton distributions. In Figs. 1 and 2 we compare the VFNS with the 

FFNS predictions for F£(x, Q 2 ) and F%(x, Q 2 ), respectively, usingTj/i 2 = Q 2 + 4m 2 b for the 

FFNS although our results are not very sensitive to this specific choice of the factorization 

and renormalization scales. As usual, fi 2 = Q 2 in the VFNS. Notice furthermore that 

the NLO- VFNS predictions for xh (short-dashed curves) are very similar to the ones for 

(2e 2 h )~ 1 F^ (dashed curves) despite the fact that (2e 2 h )~ 1 F^ = {l+a s C q )®h+\a s C g ®g, i.e. 

the 0(a s ) quark and gluon contributions almost cancel. As expected [6] the discrepancies 

Notice that here and in the following /i^ = fip = [i where /i/j and lip are the renormalization 
and factorization scales, respectively. This choice is dictated by the fact that our (and all other presently 
available) parton distributions were determined and evolved with fiR — fip, i.e. with the commonly adopted 
standard evolution equations. 
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between the predictions for xh(x,Q 2 ) in the VFNS and for (2e^) _1 F 2 h (x, Q 2 ) in the FFNS 
in the relevant kinematic region (small x, large Q 2 ) never disappear and can amount to 
as much as about 30% at very small-x, even at W 2 = Q 2 {\ — 1) far above threshold, i.e. 
W 2 ^> W 2 h = (2rrih) 2 . This is due to the fact that here the ratio of the threshold energy 
W t h = V~Sth of the massive subprocess (7*^ — > hh, etc.) and the mass of the produced 
heavy quark yfsth/rrih = 2 is not sufficiently high to exclude significant contributions from 
the threshold region. Even for the lightest heavy quark, h = c, such non-relativistic 
(/3 C = \p c \/E c < 0.9) threshold region contributions to F£(x, Q 2 ) are sizeable for W 2 £ 10 6 
GeV 2 due to significant (3 C < 0.9 contributions, and the situation becomes worse for h = b 
(cf. Fig. 4 of [6]). This is in contrast to processes where one of the produced particles is much 
heavier than the other one, like the weak CC contribution [7J |8] W + g — > tb to F% . Here 



y/s^h/m b = {m t + m h )/m b 3> 1 and the extension of the threshold region where /% ^ 0.9, 

being proportional to rrib/y/sth 1, is strongly reduced as compared to rrih/(2mh) = 0.5 
in the former case of hh production. Thus the single top production rate in W + g — > tb 
is dominated by the (beyond-threshold) relativistic region where /% > 0.9 and therefore 
is expected to be well approximated by W + b — ► t where b is an effective massless parton 
within the nucleon. In Fig. 3 we compare the LO FFNS 0E] predictions for (x,Q 2 ) 
with the corresponding VFNS ones for £6(£, Q 2 +m 2 ) where the latter refers to the W + b — > t 
transition using the slow rescaling variable [9] ^ = x(l + m 2 /Q 2 ) with m t = 175 GeV. For 
F 2 c ^(x, Q 2 ) we used fi 2 R = fi 2 F = fi 2 = Q 2 + (m t + m b ) 2 . (Notice that the fully massive NLO 
FFNS QCD corrections to W + g — > tb are unfortunately not available in the literature.) As 
expected the differences between the two schemes are here less pronounced than in the case 
of cc and bb electroproduction in Figs. 1 and 2. These results indicate that one may resort 
to the simpler VFNS with its massless h(x, fi 2 ) distributions to estimate rather reliably the 
production rates of heavy quarks, gauge bosons, Higgs scalars, etc. at Tevatron and LHC 
energies. 

As a next test of these VFNS distributions we therefore turn to hadronic W production 
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and present in Fig. 4 their NLO predictions for a(pp — > W X) as compared to the data 
PIH HI 021 H3l E] and to predictions based on the NLO CTEQ6.5 distributions [15]. Also 
shown in this figure is a comparison of our LO FFNS and VFNS predictions. Although 
quantitatively slightly different, the dominant light quark contributions in the FFNS (ud — > 
W + , us — > W + , etc.) are due to the same subprocesses as in the VFNS, but the relevant 
heavy quark contributions have been calculated via gs(d) — > cW + , gu — ► bW + , etc. as 
compared to cs(d) — > W + , bu — > W + etc. in the VFNS. Here we again expect the VFNS 
with its effective massless 'heavy' quark distributions h(x,fi 2 ) to be adequate, since non- 
relativistic contributions from the threshold region in the FFNS are suppressed due to 
y/tkh/m^b ~ M w /m C)b ^> 1. The LO gluon induced heavy quark contributions to W ± 
production in the FFNS are obtained from a straightforward calculation of the differential 
cross section [TO] da(s)/di which yields 



V2 6 s 




2 

+ (1 - 2A + 2A 2 ) In 



A 1 + 7A) (2) 



, , , 1-A + VX 



A - V A 



2 

m c 



where 



_ (m c + M w ) 2 



_ K - M wf 



a s (yU 2 ) = 47r/[9 1ii(/4 2 /Al 3 q)] and the relevant CKM matrix element(s) V qq > are taken from 
[17] . The corresponding total W hadronic production cross section relevant for Fig. 4 is 
then given by 

-1 /•! 



a 



^ cW±x (s) = I d Xl l dx 2 [g( Xl , fi 2 )s(x 2 , n 2 ) + (1 <- 2)] a(x lX2 s) (3) 

J T J t/xi 

where s(x,/x 2 ) = s(x, /z 2 ) with /i 2 = O(M^) and r = (m c + Mvy) 2 /s. Unfortunately, the 
NLO(MS) corrections to this (massive quark) FFNS cross section are again not available in 
the literature. Only quantitative LO and NLO results for the analogous process gb — > tW~ 
have been presented in [IB], but questioned in [19]. Here we just mention that we fully 
confirm the LO results for Wt production obtained in [19] at Tevatron and LHC energies. 



Taking into account that the K = NLO/LO factor is expected [19] to be in the range of 
1.2 - 1.3, our LO-FFNS predictions in Fig. 4 imply equally agreeable NLO predictions as 
the (massless quark) NLO-VFNS ones [20J shown by the solid and dashed-dotted curves 
in Fig. 4. 

In Table 1 we present our VFNS and FFNS predictions for W production at LHC and 

compare the relevant subprocess contributions to a(pp — > W ± X) at ^/s = 14 TeV. The 

light quark fusion contributions in the ud and us sector turn out to be somewhat larger 

in the FFNS than in the VFNS which is due to the fact that the u, d, s (and the gluon) 

distributions are evolved for fixed rif = 3 in the FFNS. More interesting, however, are the 

subprocesses involving heavy quarks. Here the LO-VFNS predictions are compatible, to 

within less than 15%, with the LO-FFNS predictions based on the gluon induced fixed 

order in a s subprocesses gu — » bW, gd — > cW and in particular on the sizeable CKM 

non-suppressed gs — > cW contribution. As mentioned above, the NLO corrections to 

these latter heavy quark contributions cannot be calculated for the time being. However, 

since these contributions amount to about only 15% of the total FFNS results for W ± 

production (being dominated by the light ud and us fusions in Table 1), we can safely 

employ the expected [19] K factor of K ~ 1.2 for the relevant gd — > cW and gs — > 

cW LO contributions in Table 1 for obtaining the total NLO-FFNS predictions without 

committing any significant error. The resulting total rate for W + + W~ production at LHC 

of 192.7 ± 4.7 nb is comparable to our NLO-VFNS prediction in Table 1 of 186.5 ± 4.9 nb 

where we have added the ±lcr uncertainties implied by our dynamical parton distributions 

This latter prediction reduces to 181.0 nb when using the smaller scale /i 2 = M^/4. 

The scale uncertainties of our predictions are defined by taking Mpv/2 < ji < 2Mw, 

using M w = 80.4 GeV, which gives rise to the upper limits (/i = 2M W ) and lower limits 

(n = M\y/2) of our predicted cross sections. In this way we obtain the following total 

2 Using 'standard' FFNS parton distributions [T] instead of the dynamical ones for generating the VFNS 
distributions, the dynamical NLO-VFNS prediction of 186.5 nb slightly increases to 190.7 nb. 
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uncertainty estimates of our NLO predictions at LHC: 



a(pp — > W + + W~ + X) 



186.5 ± 4.9 pdf tH (scale 



nb , 



VFNS 



(4) 



192.7±4.7 p df _ 4 ; 8 (scale 



nb , 



FFNS 



and, for completeness, at LO 



a{pp -^W + + W~ + X) 



162.1 i 3.9 p( jf ^21.8 Iscalc 



nb , 



VFNS 



(5) 



166.7 ±4.0 



I +17.3 I 
pdf —19.0 I scale 



nb 



FFNS 



where the subscript pdf refers to the la uncertainties of our parton distribution functions 
PP. For comparison, the NLO-VFNS prediction of CTEQ6.5 [15] is 202 nb with an uncer- 
tainty of 8%, taking into account a pdf uncertainty of slightly more than 2a. Similarly, 
MRST [21] predict about 194 nb. From these results we conclude that, for the time being, 
the total W production rate at LHC can be safely predicted within an uncertainty of 
about 10% irrespective of the factorization scheme. 

It is also interesting to study the dependence of the FFNS predictions for the contribu- 
tions to W production involving heavy quarks on the chosen scale \x as shown in Figs. 5 
and 6. In these figures we compare the gs — ► cW initiated production rates in the FFNS 
with the quark fusion cs — ► W ones in the VFNS and similarly the gd — > cW ones with 
the cd — > W fusion, respectively. These factorization scheme dependencies are rather mild 
for the LO-FFNS predictions, in contrast to the situation for the LO-VFNS predictions 
which stabilize, as expected, at NLO. The mild \x dependence is similar to the situation 
encountered in tW production [19J via the subprocess gb — > tW~ . 

A similar situation where the invariant mass of the produced system sizeably exceeds 
the mass of the participating heavy quarks is encountered in (heavy) Higgs H production 
accompanied by two heavy 6-quarks, for example. Here H = Hg M ; h°, H°, A denotes 
the SM Higgs-boson or a light scalar h°, a heavy scalar H° and a pseudoscalar A of 
supersymmetric theories with Mh ^100 GeV. In the FFNS the dominant production 
mechanism starts with the LO subprocess gg — > bbH (qq — > bbH is marginal), to be 
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compared with the bb fusion subprocess in the VFNS starting with bb — > H at LO. Again, 
y/sth/iTib — (2m;, + Mff)/m(, ~ Mnjm^ ^> 1 in the FFNS which indicates that the simpler 
LO and NLO(NNLO) VFNS bb fusion subprocesses do provide reliable predictions. Within 
the scale uncertainties it turns out that the FFNS and VFNS predictions at NLO are 
compatible [221 1231 121 125!, usin g the MRST2002 and CTEQ6 parametrizations of the 
relevant parton distributions [2T| l25]. This result holds for scale choices \i^p = (| to ^)\/%h 
with \/Jth/^ being considered as a suitable 'central' choice in b(x, fi 2 F ) for calculations based 
on the bb fusion process in the VFNSj^] It should, however, be mentioned that the VFNS 
rates somewhat exceed [22} |23} EH [25] the corresponding FFNS Higgs-boson production 
rates by about 3 10-20%. 

Finally let us note that all our results and comparisons concerning the VFNS hold 
irrespective of the specific parametrizations used for the 'heavy' h(x, fi 2 ) distributions: 
when comparing our VFNS distributions, generated from using our dynamical distribu- 
tions [1] as input, with the ones of CTEQ6 [26J or CTEQ6.5 [15] the relevant ratios 
c(x,M^) CTE q/c(x,M^) gjr _vfns and b(x, m 2 ) C TEQ/b(x, m 2 ) Gm _ YFm vary, for 10 -4 < 
x 0.1, at most between 0.9 - 1.1 at LO and NLO. Similar results hold when using 
other VFNS distributions, e.g., those of [Sj. This is illustrated more quantitatively in 
Fig. 7 where we compare our c- and 6-distributions, together with the important gluon- 
distribution, with the ones of CTEQ6 [26] and CTEQ6.5 [15] in the sea- and gluon-relevant 
x-region x £ 0.3 at Q 2 = M^. The ratios for the light u- and ^-distributions are even 
closer to 1 than the ones shown in Fig. 7, typically between 0.95 and 1.05 which holds 
in particular for the CTEQ6 distributions when compared to our ones. Incidentally the 
VFNS under consideration and commonly used J5j |26] is also referred to as the zero-mass 
VFNS. Sometimes one uses an improvement on this, now known as the general-mass VFNS 
[151 1211 [23 123 [291 130] , where mass-dependent corrections are maintained in the hard cross 
sections. This latter factorization scheme interpolates between the strict zero-mass VFNS, 

3 Thc independent variation of \i,f and considered in |221 1231 1241 125] is, as mentioned before, not 
strictly compatible with the utilized parton distributions determined and evolved according to \ir = fip . 
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used in our evolution to Q 2 ^> m\, and the (experimentally required) FFNS used for our 
input at Q 2 = m\. As expected and shown in Fig. 7, scheme (input) differences at lower 
Q 2 = O(m^) only marginally affect the asymptotic results at Q 2 = 3> m 2 b where 
the CTEQ6.5 parametrizations [f 5J (corresponding to a general-mass VFNS) become very 
similar to the ones of CTEQ6 [26] and our GJR-VFNS (corresponding to the zero-mass 
VFNS). As stated repeatedly before, this is essentially due to the dominance of the large 
evolution effects over the minor differences involved at the lower scales, e.g. at Q 2 = 0(m\). 
These asymptotic similarities are particularly relevant for the simplified (vanishing m C) &) 
calculations of the production rates of very massive particles where massive c- and 6-quark 
threshold region contributions are strongly suppressed. 

To summarize, we generated radiatively two sets of VFNS parton distributions, based 
on our recent LO and NLO dynamical parton distributions pQ obtained in the FFNS. 
Within the VFNS additional heavy quark distributions h(x, Q 2 ) = h(x, Q 2 ) are generated 
perturbatively via the common massless Q 2 -evolution equations by imposing the boundary 
conditions h(x, Q 2 = ml) = for h — c,b, t. We have confronted the VFNS and FFNS 
predictions in situations where the invariant mass of the produced system (hh, tb, cW, tW, 
Higgs-bosons, etc.) does not exceed or exceeds by far the mass of the participating heavy 
flavor. In the former case (e.g. i 7 ^ in deep inelastic cc production where \fsth/rn c = 2) the 
VFNS predictions deviate from the FFNS ones by up to about 30% even at Q 2 ;§> m 2 . 
In the latter case (e.g. F^9 in deep inelastic weak charged current tb production where 
\/Ith/rrib — mt/rrib ^> 1 these deviations are appreciably less, within about 10%, which is 
within the margins of renormalization and factorization scale uncertainties. As a further 
example of the agreement between the VFNS and FFNS predictions in situations where the 
invariant mass of the produced system far exceeds m c ^ we considered the corresponding 
W^- boson production rates at the Tevatron and at the large hadron collider (where e.g. 
y/s t h/m c ,b — Mw/m c ,b ^> 1 for cW and bW production, respectively). For y/s = 14 TeV 
the NLO-FFNS predicts a(pp -> W + + W~ +X) ~ 192.7 nb with an uncertainty of 5%, to 
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be compared with the NLO-VFNS prediction of 186.5 nb and an uncertainty of 6%. The 
cited uncertainties include also the scale uncertainties due to varying the renormalization 
and factorization scales //r = //f between and 2M^. (It should be emphasized again 

that the scale choice /z# = \lf is dictated by all presently available parton distributions 
which have been determined and evolved according to /z# = //p.) Furthermore, a similar 
agreement is obtained for hadronic (heavy) Higgs-boson production when the dominant 
FFNS subprocess gg — ► bbH (where \/§th/mb = (2rrib + M H )/mb ^> 1) is compared with 
the VFNS 6-quark fusion subprocess (66 — > H, etc.). 

These results indicate that the simpler VFNS with its effective treatment of heavy 
quarks (c, 6, t) as massless partons can be employed for calculating processes where the 
invariant mass of the produced system is sizeably larger than the mass of the participating 
heavy quark flavor. The uncertainty of such calculations is process (and somewhat energy) 
dependent when compared with the predictions of the FFNS where the effects of finite 
heavy quark masses are nowhere neglected. Taking into account the uncertainties of parton 
distributions and scale choices as well, the total W production rate at LHC can be 
predicted within an uncertainty of about 10%, which becomes significantly smaller at the 
Tevatron. Similarly the Higgs production rates at LHC are predicted with an uncertainty 
of 10-20% where the VFNS production rates exceed the FFNS ones by about 20% at LHC. 

A FORTRAN code (grid) containing our new LO and NLO(MS) light (u,d,s;g) 
and heavy (c, 6, t) parton distributions in the VFNS, generated from our recent 
dynamical ones in the FFNS [I], can be obtained on request or directly from 
http : // doom . physik . uni-dortmund . de . 

Acknowledgements 

We thank J. Campbell for a clarifying correspondence. This work has been supported in 
part by the 'Bundesministerium fiir Bildung und Forschung', Berlin/Bonn. 



9 



References 

[1] M. Gluck, P. Jimenez-Delgado, E. Reya, Eur. Phys. J. C53 (2008) 355 

[2] M. Gluck, E. Reya, A. Vogt, Eur. Phys. J. C5 (1998) 461 

[3] T. Stelzer, Z. Sullivan, S. Willenbrock, Phys. Rev. D56 (1997) 5919 

[4] J.C. Collins, W.K. Tung, Nucl. Phys. B278 (1986) 934 

[5] S.I. Alekhin, Phys. Rev. D68 (2003) 014002 

[6] M. Gluck, E. Reya, M. Stratmann, Nucl. Phys. B422 (1994) 37 

[7] M. Gluck, R.M. Godbole, E. Reya, Z. Phys. C38 (1988), 441; C39 (1988) 590(E) 

[8] U. Baur, J.J. van der Bij, Nucl. Phys. B304 (1988) 451 

[9] R.M. Barnett, Phys. Rev. Lett. 36 (1976) 1163; 

R. Barbieri, J. Ellis, M.K. Gaillard, G.G. Ross, Nucl. Phys. B117 (1976) 50 

[10] C. Albajar et al, UA1 Collaboration, Z. Phys. C44 (1989) 15 

[11] J. Alitti et al., UA2 Collaboration, Phys. Lett. B276 (1992) 365 

[12] B. Abbott et al, D0 Collaboration Phys. Rev. D61 (2000) 072001 

[13] F. Abe et al, CDF Collaboration, Phys. Rev. Lett. 69 (1992) 28 

[14] A. Abulencia et al., CDF Collaboration, J. Phys. G34 (2007) 2457 

[15] W.K. Tung et al., CTEQ Collaboration, JEEP 02(2007)053 

[16] F. Halzen, C.S. Kim, Int. J. Mod. Phys. A2 (1987) 1069 

[17] W.-M. Yao et al., Particle Data Group, J. Phys. G33 (2006) 1 

[18] S. Zhu, Phys. Lett. B524 (2002) 283; B537 (2002) 351(E) 

10 



[19 
[20 
[21 
[22 
[23 
[24 

[25 
[26 
[27 
[28 



J. Campbell, F. Tramontano, Nucl. Phys. B726 (2005) 109 

R. Hamberg, W.L. van Neerven, T. Matsuura, Nucl. Phys. B359 (1991) 343 

A.D. Martin et al., Eur. Phys. J. C23 (2002) 73; Phys. Lett. B531 (2002) 216 

S. Dittmaier, M. Kramer, M. Spira, Phys. Rev. D70 (2004) 074010 

S. Dawson, C.B. Jackson, L. Reina, D. Wackeroth, Phys. Rev. Lett. 94 (2005) 031802 

J. Campbell et al., 'Physics at TeV Colliders' (3 rd Les Houches Workshop, 2003), 



|hep~ph/0405302l 

S. Dawson, C.B. Jackson, L. Reina, D. Wackeroth, Mod. Phys. Lett. A21 (2006) 89 
J. Pumplin et al., CTEQ Collaboration, JEEP 07 (2002) 012 
S. Kretzer et al., Phys. Rev. D69 (2004) 114005 
M.A.G. Aivazis et al., Phys. Rev. D50 (1994) 3102 



[29] M. Buza et al., Eur. Phys. J. CI (1998) 301; 

A. Chuvakin, J. Smith, and W.L. van Neerven, Phys. Rev. D61 (2000) 096004 

[30] R.S. Thorne and R.G. Roberts, Phys. Rev. D57 (1998) 6871; 
R.S. Thorne, Phys. Rev. D73 (2006) 054019 



% 



\ 

b 



O 
fa 

o 
fa 

2 

co 
Z 
fa 
fa 



O 
fa^ 

o 
fa 

z 

co 
Z 
fa 
> 



i 

+ 

+ 



+ 



3 



X 

13 



fa 

T 

pi 

bJ3 



o 

T 



o 

T 



t3 



CD 

fa 
CO 
i — I 

CS 
fa 



cd 
fa 

CD 



OS 

i — I 

GO 

CO 

CO 
OS 



00 
UO 



o 

1 — I 

X 

CO" 

co 



o 

CO 



CO 
1 — I 
CN 



CO 



o 

1 — I 

X 



O 
CM 



CD 

o 



CO 
CM 

fa 



o 

1 — I 

X 
CD 



CD 
O 



(130.3 


(4.5) 


2) x 10- 


i — 1 

co 


(24.3) 


CO 


CO 


os 


CD 




06 
i — i 


to 


CD 

as 


CO 


OS 

CM 






-f 






(52.6) 


CN 
i — 1 


2)xl0" 


(2.0) 


(12.2) 


CD 


CO 


CN 


CO 




O 
CD 


1 — 1 


CO 


CM 


i — i 



o 

1 — I 

X 

CM 

OS 
CM 



o 

1 — I 

X 



CD" 




o 

1 — 1 




cm" 


fa 




X 


i — I 


CM 




CO 


H 


i — 1 






o" 








CS 


fa 


CO 


l> 


fa 


CO 




i — 1 


fa 


00 




CO 




i — 1 






fa 






T3 


cc 


fa 




CO 


pi 




13 


o 


o 



o 

1 — I 

X 

1 — I 

UO 



fa 

o 



O 



CD 
CD 



cm 

05 



CD 

OS* 
CD 



CM 

i — I 

00 

l\ 



OS 



CM 
CD 



CD 



Co 

fa 

CD 



OS 



os 



o 



o 

+3 



CP CO 

a3 o 

CO 

PI pi 

■2 2 

Pi +^ 

fa .2 

■C ^ 



3^ 



CJ 

o 

ft 

o 



to fa 
fa f- 1 

ci3 

§ & 



53 o3 

^ .g 

ft ai 
cp 
fa 

- 

cp 



co 
Z 
fa 
fa 

o 
fa 
H 

d 
fa 



pi 

.2 

+j 
o 
pi 

o 

ft 

-H 
cc 

PJ 

O 



o 

fa 

bJ3 

a 

xs 

_p| 

co 
Z 
fa 
> 
o 
fa 



ft CO 

O fa 

a3 

CJ 1 



x 
o 
fa 

'3 
+j 
sh 

o 
o 
P! 
Pi 

o 
fa 
H 



H XJ 

.S aj 

J CJ 

8 2 

0) " 

fa CD 

fa 

CO > 

CP r-H 

+= O 

ce > 



8 ^co 

S ^ £ 
° fa 
fa 



a 
fa 

bJO 
PI 



bJO 
PI 



O CP 

PI . . 

bC ^ 

o 

CO "+h 

fa 
fa 



x 
o 
+= 

o 
fa 



T3 
o 

ft 



T3 
Pi 

ct 

CO 

Z 
fa 

> 



CO 

'cc 

-. 

fa 

o 

fa 



fa 

CP 

O fa « 



^ p 



Pi 

a3 



fa 



3. 



o 
fa 

z ^ 

^ CP 
I — I 

Co t^j 

O fa 
^ PI 
CP -rH 

fa PI 

H I 

^fa 'bb 

pi 

■x) P3 

CO ^ 

si 

CO 

CO CP 

C fa 

CP 

•P £ 

CM 
CP . 
N t — I 



fa ^ H <+H 

z a , , I c 



o 11 

H ^ ^ 

» s a. 



o 
2 2 

CO ^ 
CP , 



3.5 



- 1 1 — l l l l l 1 1 1 1 — l l l l l 1 1 1 1 — l l l l l l 



- 1 1 — l l l l 1 1 1 1 1 — l l l l 1 1 1 1 1 — I — I I I 1 1 

XC — : 



8 h 2 _ 



— FFNS 

- - VFNS 



\ W = Q 2 (GeV 2 ) 




NLO 




lO'VlO" 4 



Figure 1: The predicted x-dependencies of ^F£(x, Q 2 ) in the FFNS and VFNS at some 
typical fixed values of Q 2 . For the FFNS the renormalization and factorization scales are 
chosen to be [i 2 R = [i 2 F = fi 2 = Q 2 + Am 2 with m c = 1.3 GeV, and, as usual, fi 2 = Q 2 
for the VFNS. The NLO-VFNS charm distribution is given by xc(x, Q 2 ) as shown by the 
short-dashed curves. 




Figure 2: As in Fig. 1 but for bottom production, i.e. |F 2 b (x,Q 2 ), choosing fi 2 R = \j? F = 
/i 2 = Q 2 + 4m 2 with m fe = 4.2 GeV for the FFNS. The short-dashed curves show the 
NLO-VFNS bottom distribution xb(x,Q 2 ). 




Figure 3: LO predictions for the ^-dependencies of the weak charged current structure 
function \F^2( X ->Q 2 ) f° r tb production in the FFNS at some typical fixed values of Q 2 . 
The momentum scale is chosen to be /jl 2 r — n 2 F = pi 2 — Q 2 + (m t + m b ) 2 with m t = 175 
GeV. These predictions are compared with the bottom distribution £6(£, Q 2 + m 2 ) in the 
VFNS where £ = x(l + m 2 /Q 2 ). 
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Figure 4: Predictions for the total W + + 14 7- production rates at pp colliders with the 
data taken from PH Ell H21 EH E]- The LO and NLO GJR parton distributions in 
the VFNS have been generated from the FFNS ones [1] as described in the text. The 
NLO-VFNS CTEQ6.5 distributions are taken from [15]. The adopted momentum scale 
is fiR = [If = A* — Mw- The scale uncertainty of our NLO GJR predictions, due to 
varying /z according to \Mw < ^ < 2Mw, amounts to less than 2% at yfs = 1.96 TeV, 
for example. The shaded region around our central GJR predictions is due to the ±lcr 
uncertainty implied by our dynamical NLO parton distributions [T]. 
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Figure 5: The scale dependence (/xr = fip = fi) of the LO-FFNS contribution to the total 
W + + W~ production rate due to the subprocess gs — > cW compared to the LO and NLO 
ones in the VFNS due to cs — > W fusion. The results refer to the pp-LHC (^/s = 14 TeV) 
and to the pp-Tevatron (y/s = 1.96 TeV) with the latter ones being multiplied by a factor 
of 10 as indicated. 
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Figure 6: As in Fig. 5 but for the FFNS subprocess gd — > cW to be compared with cd W 
in the VFNS. The results for the Tevatron (^/s = 1.96 TeV) are multiplied by a factor of 
5 as indicated. 




Figure 7: Comparing our present (GJR-VFNS) dynamical parton distributions generated 
in the VFNS at NLO(MS) with the ones of CTEQ6 [25] and CTEQ6.5 [B] at Q 2 = M^. 



